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We analytically and numerically investigate the ground state of the spin-orbit coupled spin-1 Bose- 
Einstein condensates in an external parabolic potential. When the spin-orbit coupling strength 
n is comparable with that of the trapping potential, the density distribution centers of different 
components of the spinor condensate deviate evidently from the trap center in the plane wave and 
stripe phases. When k ^> 1, the magnitude of this deviation decreases as k is getting larger and 
larger. Correspondingly, periphery half-skyrmions textures arise. This deviation can be reflected by 
the non-uniform magnetic moment in the z direction, T z . With the manipulation of the external 
trap, the local magnitude of T z can be increased evidently. This kind of increase of T z is also 
observed in the square vortex lattice phase of the condensate. 
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I. INTRODUCTION 

The spin-orbit (SO) coupling is an essential mecha- 
nism for most spintronics devices, and leads to many fun- 
damental phenomena in condensed matter physics and 
atomic physics. For example, SO coupling gives rise to 
quantum spin Hall effect in electronic condensed matter 
systems [l|. Recently, artificial external Abelian or non- 
Abelian gauge potentials coupled to neutral atoms have 
been generated by controlling atom-light interaction 
[||, which provides the possibility of spin control in neu- 
tral atom systems through laser fields. Without the SO 
coupling, spinor Bose-Einstein condensates (BECs) have 
received extensive explorations, including the spin dy- 
namics and textures [6-20]. In the presence of the SO 
coupling, the neutral atoms in the effective abelian and 
non-abelian gauge fields behave like electrons in an elec- 
tro magnetic field or electrons and abundant phenom- 
ena arise. Thus, SO coupled cold atom systems have 
attracted intensive attention since last few years (2ll-[39j . 
For example, the condensate may become self- trapped 
as a result of the spin-orbit coupling and the nonlin- 
earity, resembling the so-called chiral confinement (33| . 
Dynamical oscillations of the condensate with SO cou- 
pling are studied in [34j], where oscillation period, sim- 
ilar to the Zitterbewegung oscillation, was found. Be- 
sides, SO coupled BECs with dipole-dipole interactions 
[36l] an d rotating trap [37|, [38[, have also been studied. 
In [39[, SO coupled atomic spin-2 BEC was investigated, 
and square or triangular density patterns were reported. 
Both patterns evolve continuously into striped forms with 
increased asymmetry of the SO coupling. 

The SO coupled spinor condensate develops a spon- 
taneous plane wave phase or stripe phase because of 
the interaction energy [40l-[42|. In this case, the har- 
monic length of the external trap is much greater than 



the wavelength of the stripe, indicating a weak external 
trap compared to the SO coupling strength. In two com- 
ponents condensate, the phase diagram including half- 
quantum vortex, vortex lattice and stripe states has been 
obtained in Ref. [23j . For the case with strong external 
traps and strong spin-orbit coupling, half-quantum vor- 
tex state and Skyrmion lattice patterns can arise [24j , 
where the atomic interaction is relatively weak. How- 
ever, spin-1 BECs in the presence of external trap with 
neither too strong nor too weak SO coupling strength, 
i.e., k ~ 1, are less studied. The atom-atom interac- 
tion needs to be large enough to produce trapped plane 
wave or stripe state in an external trap. In the trapped 
plane wave or stripe state, there exist nontrivial detailed 
structures, as we will see in the following sections. 

In the present paper, we both analytically and numeri- 
cally investigate certain characteristics of the spin-1 BEC 
with SO coupling strength n ~ 1 in external parabolic po- 
tential. Comparing with the homogenous case, we find 
that the ground state wave functions of different compo- 
nents of the condensate dislocate and exhibit finite sep- 
aration, which is reflected by the non-uniform magnetic 
moment distribution in z direction. In particular, we 
find associated hidden spin textures in the usual plane 
wave and stripe phases. Also, effects of the manipula- 
tion of external potentials on the magnetic moment in z 
direction, are explored in the trapped plane wave, stripe 
and the square lattice structure phases. The separation 
in the trapped stripe state is more complex. By apply- 
ing an anisotropic external potential, the trapped stripe 
state develops a spin texture characterized by aligned 
half-skyrmions. 

The rest of the paper is organized as follows. In Sec. 
HH we propose an variational approach to study the sep- 
aration of density distributions between different compo- 
nents in the plane wave phase. In Sec. we numer- 
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FIG. 1: Schematic picture of the band structure with SO 
coupled cold atoms. 



ically confirm the analytical results obtained in Sec. ITT1 
and exhibit the hidden spin texture in the plane wave 
phase. Sec. [TV] focuses on the hidden spin textures in 
stripe phase. Finally, the conclusions are summarized in 
Sec. ED 



II. THEORETICAL FRAMEWORK 

We consider the quasi-2 dimensional spin-1 BECs with 
Rashba SO coupling in the x — y plane confined in an 
external parabolic potential. Without regarding to the 
interaction, the single particle part of the Hamiltonian 
without external potential is 



dr 



2/cp • F 



2M 



(1) 



where F = (F x , F y ) is spin-1 representation of Pauli ma- 
trices. The SO coupling strength n is experimentally re- 
lated to the wavelength of the laser beams. The three 
components spinor is \l> = (^i , ^o, ^-i) T (the super- 
script T stands for the transpose) and M is the atomic 
mass. 

In order to gain some intuition, we first focus on the 
homogeneous condensate. Thus, the single-particle spec- 
trum of the Hamiltonian (pQ) has the following three 
branches with different helicity in the momentum space 
(see Fig. [Q). 
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The single-particle ground state is in the negative helicity 
branch with po = |p| = y/2K, and the wave function is 
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FIG. 2: (Color online) The SO-coupling-energy function H s 
of (a) Si and (b) S-i with r\ = 0.0096 and k = 2.0. The 
light-red balls at points 5i = S-i =0 feel forces pointing to 
the equilibrium position. 



given by 
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where </> p = a.Tg(p x + ip y ). All the states with the same 
|p| are circularly degenerate. 

Taking account of the interaction and the external po- 
tential, the mean-field version of the total Hamiltonian 
is 
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where the coupling constants cq and characterize the 
density-density and spin-exchange interactions, respec- 
tively. The length, time and energy are respectively 
scaled in units of ^/h/Muo, 1/uj and hw. a; is the radial 
trapping frequency. 

To investigate this hidden separation of the plane wave 
phase in the direction perpendicular to the propagation 
of the plane wave, we suppose a plane wave propagating 
along x axis and use the ansatz: 
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where r] and £ are positive dimensionless variational pa- 
rameters related to the width of the condensate, and 
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FIG. 3: (Color online) (a) Dependence of 6 on 77 with k — 
2.0. (b) Dependence of 5 on k with 77 = 0.0096. The black 
dotted line corresponds to the case when only the SO part 
of the mean field energy is minimized. The solid red line 
corresponds to the case of counting in both the SO and the 
external potential parts of the mean field energy. 
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- is the normalization coefficient. Here 

we introduce the deviation S m from the origin in the y 
direction. A deviation in the x direction can also be sup- 
posed and this deviation is found to be zero to minimize 
the total mean field energy. 

Without the SO coupling term in the Hamiltonian, the 
condensate feels comfortable when the distribution cen- 
ters of the three components occupy the center of the 
external parabolic potential. It is not the same situa- 
tion any longer once the SO coupling is included. In the 
presence of the SO coupling, the gaussian-profile density 
distribution with Si = So = S-i = is not the ground 
state . We can introduce the effective virtual force to de- 
termine the deviating direction of gaussian-profile den- 
sity of each component. Thus, the mean- field energy 
would be minimized. We define the effective virtual force: 
fi = —dH s /dSi\s i= o (i = 0, ±1). After substituting the 
ansatz in Eq. (j3j) into the SO part of the total energy, 
we have 
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Since 77, £ and n are positive parameters, the direction 
of the dragging force have been fixed. The component 
with magnetic quantum number m = 1 feels a dragging 
force pointing to the left side of the propagating plane 
wave direction and the component with magnetic quan- 
tum number m = — 1 feels the same force but in opposite 
direction, while no force acts on the components with 
magnetic quantum number m = 0. Based on these effec- 
tive virtual forces, we suppose Si = — S-i = S, So = 0, 
and £ = 77. In fact, further calculations show that the 
deviation is not sensitive to the relative width of the con- 
densate for the present concerned range of values of the 
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FIG. 4: (Color online) (a) The density (upper panels) and 
phase (bottom panels) distributions of the plane wave phase 
with k = 2.0, Co = 500.0 and C2 = —5.0. (b) The non-uniform 
density of magnetic moment along the z direction, T Z: in the 
plane wave phase. 



SO coupling strength. After substituting the ansatz in 
Eq. ^ into the SO part of Eq. (|2]), we have 



-^e-^ 2r]Po2d2 V2(po + VPo 2 S) 



(5) 



Starting from H S1 we can find that the minimum func- 
tional energy requires the condition: 



-1 + Vl + 4r/ 



>0. 
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In Fig. [2j the SO part of the total mean-field energy with 
respect to Si (Si = S) and 6-1 (S-i = —S) is plotted. The 
light-red balls, which deviate from the equilibrium posi- 
tion feel forces as expressed in Eq. (jl]). We know that 
spin currents accompany charge currents in high-mobility 
two-dimensional electron systems (43| , and thus spin ac- 
cumulation will be induced with the presence of the SO 
coupling. Starting from the single particle Hamiltonian 
in Eq. (pQ) , we can obtain the force felt by the components 
of the condensate in the Heisenberg picture by following 
calculations in Ref. [44j |: 



f = 2K 2 /(HM 2 )(pxe z )F z . 



(7) 





FIG. 5: (Color online) (a) As the condensate is being 
squeezed, the peak value of the T z increases with respect to 
(jJ±/(jJ, where uj± is the trapping frequency in the squeezing 
direction, (b) The corresponding spin texture, showing the 
hidden half-skyrmions in the periphery of the condensate in 
the plane wave phase. The relevant parameters are k = 2.0, 
uj ± /uj = 4.0, c = 500.0 and c 2 = 5.0. 



Here we should mention that, the above introduced effec- 
tive virtual forces in Eq. @ come from small deviations 
from the ground state (equilibrium position), which is 
not the same as the true force of the particles during the 
dynamical motions, as expressed in Eq. ([7j). 

The dependence of S on r] and n is plotted in Fig. [3] 
From Fig. [3fa) with n = 2.0, we can find that the de- 
viation does not change dramatically with respect to r] 
and chooses values 0.2 ~ 0.35. More calculations show 
that the deviation function of SO coupling strength is 
5 = l/y/2K as r] — >• (meaning an unlimited condensate 
width corresponding to the case of zero external poten- 
tial). This limit behavior is not affected by the interac- 
tions between particles. Figure [3^b) shows the depen- 
dence of S on k with r] = 0.0096. For the black dotted 
line, only the SO part of the mean field energy is mini- 
mized. As k increases, i.e., is approaching the strong SO 
coupling regime, the deviation diminishes at a rate pro- 
portional to From Fig. [3f b), we can see from the 
black dotted line that the deviation diverges as the SO 
coupling strength getting weaker and weaker ( k — >• 0), 
which is nonphysical. This divergence can be canceled 
by counting in the external potential, as shown by the 
solid red line of Fig. [3fb). 

To minimize the mean-field energy, there needs to be 
deviation of density distribution between different com- 



ponents of the condensate, resulting in polarized mag- 
netic momentum distribution {Tz). It is necessary to no- 
tice that the polarized magnetic momentum distribution 
has nothing to do with spin-dependent force expressed in 
Eq. ([7j). Through minimizing the mean- field energy, the 
polarized magnetic momentum distribution in the plane 
wave phase is the consequence of the cooperation of all 
the components of the condensate in the presence of the 
SO coupling and the external trap. 

For the stripe phase of the condensate, the matter wave 
is composed of two coherent wave vector states that prop- 
agating face to face along the x axis. We can choose the 
ansatz: 
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(8) 

where C = po is the normalization coefficient. Start- 
ing from the SO part of Eq. ((2]), the deviation formula 
similar to Eq. (|6j) can be found: 



0. 



(9) 



Equation (|9]) requires 5 = 0, i.e., the deviation which oc- 
curs in the plane wave phase do not reside in the stripe 
phase. However, as we point out in section IIV[ there is 
actually small deviation in the stripe phase due to the 
atomic interaction. As the external potential becomes 
more and more anisotropic, this deviation manifests it- 
self more clearly. Deviations in anisotropic external po- 
tential is more complex for the stripe phase and we need 
to consider an additional degree of freedom to describe 
the deviation, as we will introduce in Sec. HVl 
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III. NUMERICAL RESULTS IN THE PLANE 
WAVE PHASE AND EXPLORATIONS BY 
MANIPULATION OF THE EXTERNAL 
POTENTIAL 

In this section we numerically investigate the devia- 
tion pointed out in section [TTJ As shown in the previous 
analytical calculations, the deviation S is small and it is 
usually difficult to identify the existence of this kind of 
deviation by only keeping eyes on the density distribu- 
tion. In Fig. HJa), we show the plane wave phase of 
spin-1 condensate. Theoretical explorations show that 
plane wave phase can be realized in the 87 Rb condensate 
because the spin-exchange interaction coefficient c<i < 

0. 

By resorting to the density of magnetic moment along 
the z direction: 

7- z = |Vi| 2 -|V-i| 2 , (10) 

we can specify this kind of deviation. If density distribu- 
tion of particles in component ?pi coincides with that in 
component a uniform T z would be expected. Oth- 
erwise, non-uniform T z would appear, indicating the de- 
viation of particle density distributions in components ipi 
and ip-i. We show the distribution of T z in Fig. H(b). 
From Fig. ffl[b), it is clearly shown that T z chooses pos- 
itive values in the bottom right corner of the panel, i.e., 
the density distribution of the ipi component has been 
dragged along direction (x — y)/y/2 (x and y are unit 
vectors along axis x and y, respectively). The phase dis- 
tribution in Fig. HJa) shows the plane wave propagates 
along the direction {—x — y)j\f2. The numerical results 
agree with the analytical calculations in section ITT1 

In order to compare the analytical and numerical re- 
sults, we start calculations from Eq. (|2]) with cq = 500.0 
and C2 = —5.0. From Fig. H^a), we find the deviation of 
the density peaks between components with m = 1 and 
m = —1 is 0.25. By referring to minimum of the total 
energy, we can find that such a deviation corresponds to 
rj = 0.006, i.e., a width of the condensate A = 3.22, which 
is compatible with the size of the numerically obtained 
condensate. However, as the interaction grows stronger, 
the consistency between analytical and numerical results 
gets worse and worse because the gaussian configuration 
ansatz is no longer appropriate. In fact, the density of 
each component shows asymmetrical distribution refer- 
ring to its peak density point. 

As the small deviation can be reflected by the polarized 
magnetic moment distribution, possible ways to magnify 
the local magnetic moment would be useful. We im- 
plement an anisotropic external trapping potential, and 
show the response of the local magnetic moment to this 
anisotropism. In Fig. [5ja), we show the variation of the 
peak value of T z with respect to the uj±/uj, where u± 
is the trapping frequency in the squeezing direction. As 
uj±/uj increases, the local magnetic moment gets larger 
and larger, which is helpful to assure the small deviation 
in the experiments. The spin texture, i.e. half-skyrmions, 
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FIG. 7: (Color online) The panels in the bottom are the spin 
texture obtained by squeezing the condensate in the stripe 
phase with k — 2.0, Co — 500.0, C2 — 5.0 and uo±/u — 8, 
where uo± is the trapping frequency in the squeezing direction. 
Spin texture on the right is the magnification of the squared 
part in the left. The panel in the top shows the corresponding 
phase distribution of the ipi component in the stripe phase. 

is plotted in Fig. [5(b). For the square lattice structure 
of the condensate, half-skyrmions lattices would appear, 
the squeezing of the external potential takes up the role 
of magnifying the local z component of the magnetic mo- 
ment, which is similar to that in the plane wave phase. 
We should mention that, as gets larger, the con- 

densate would cross from the square lattice phase to the 
plane wave phase. 



IV. HIDDEN SPIN TEXTURES IN THE STRIPE 
PHASE 

Through the analytical and numerical calculations, we 
know that there is deviation between the density distri- 
butions of components and ip-i for the plane wave 
phase. This deviation occurs in the direction perpendic- 
ular to the propagating direction of the plane wave. In 
this section we focus on the stripe phases. The stripe 
phase is predicted to exist in the spinor condensate of 
23 Na with the spin-exchange interaction C2 > [411 ]. 

In the following, we take the parameters n = 2.0, 
Co = 500.0 and C2 = 5.0, and the trapped stripe phase is 
found numerically. The density and phase distributions 
of component ipi are plotted in Fig. [6ja). Figure Efb) 
shows the spin texture of the stripe phase. The hidden 
half-skyrmions can be clearly seen in the periphery of the 
condensate. The spin texture of a pair of half-skyrmions 
is shown in the bottom of Fig. EJb). By applying ex- 
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ternal manipulation on the condensate, i.e., squeezing 
the condensate in certain direction, the spin texture can 
manifest themselves more clearly, which is shown by the 
panels in the bottom of Fig. The color on the arrows 
represents the local magnetic moment in the z direction 
(the modulus of the magnetic moment have been normal- 
ized to unity). As the color changes from blue to dull-red, 
the relative atom number changes between components 
tpi and tjj-i. A line of paired half-skyrmions can be seen. 
The panel on the right is the magnification of the squared 
area in the panel on the left. From the spin texture in 
Fig. [TJ there seems to be small deviation between com- 
ponents -01 and ijj-i because of the evident color contrast 
in the paired half-skyrmions line. 

In order to understand the numerical results in Fig. 
qualitatively, we introduce the following ansatz: 

( 1 _ l^^_ e ^(Fg-P0 2 ^ 2 +^(y-yo+^) 2 ) _|_ e ~W(P e -P0 2 r{x 2 +R(y-yo-$) 2 ) 

(11) ' 

where, C = po ^J^. is the normalization coefficient, and 
R can fix the ratio of the size scales of the condensate 
in x and y directions. Repeating the same calculations 
as in Sec. [Ill we can obtain the dependence of the en- 
ergy on parameters S and yo. Here S represents separa- 
tion between the two gaussian packets bearing counter- 
propagating plane waves and yo characterizes the average 
separation of packets in m = 1 and m = — 1 components. 

From Fig. Ufa), where the parabolic potential and SO 
parts of the energy is plotted as a function of S and yo , we 
can see that the energetic minimum is shifted in the neg- 
ative direction of the S axis (marked by light-red ball), 
while there is no shift in the yo axis. By counting in 
the non-zero S in the ansatz (fTTj) . the phase distribution 
shows the same configuration with that of the numeri- 
cal result in Fig. In Ref. [45], a variational ansatz 
representing a modified striped state was proposed. The 
phase distribution in the present paper is consistent with 
the particle current found in Ref. [45|. 

Thus far, the roles of the interaction have not been in- 
cluded. In Fig. [8]^b), we plot the interaction energy E int 
as a function of S and yo. There are eight energetic min- 
imums (marked by light-red ball in the quasi triangular 
contour) around the center. Therefore, the role of the in- 
teraction is to drag the energetic minimum off the center 
(5 = 2/0 = 0). The color distribution (blue to dull-red) in 
Fig. [3(b), which represents non-uniform magnetic mo- 
mentum component T z , indicates a non-zero yo has been 
chosen. 

In the end, we discuss and exhibit some key experimen- 
tal parameters. The rea l SO coupling strength in the 
present paper is 7 = K\/hmuj. A Rashba SO coupling 
can be induced in spinor 87 Rb cold atoms Q. Experi- 
mentally, a tight harmonic confinement along z is neces- 
sary to realize the pancake-like condensate. If a spinor 




FIG. 8: (Color online) The dependence of (a) the Htrap + H s 
and (b) Hint on 5 and yo, with parameters k = 2.0, r\ = 
0.0096, co = 500.0, c 2 = 5.0, and R=64. 

condensate composed of N = 10 4 atoms is realized and 
the radial and axial trapping frequencies are chosen as 
2tt x S50Hz and 2tt x 20Hz, respectively, the correspond- 
ing s-scattering lengths are ao = 68. 76a b (total spin 
channel F tota i = 0) and CL2 = 70.87clb (total spin channel 
Ftotai = 2) for 87 Rb cold atoms and ao = 138.28ob and 
a2 = 134. 21a# for 23 Na cold atoms. 



V. CONCLUSIONS 

In summary, we have analytically and numerically 
studied deviations between the density distributions of 
the ipi and ip-i components in SO coupled spin-1 BEC 
with external potential, which results in the hidden spin 
textures in the trapped plane wave and stripe phases. 
The polarized magnetic moment distribution in z direc- 
tion, T z , can mark the existence of the deviation. Manip- 
ulation of the external potential can be utilized to mag- 
nify the local T z in both the plane wave and the square 
lattice phases. Although T z in the stripe phase is small, 
nontrivial spin texture dwells in the condensate. As the 
squeezing of the external trap gets harder and harder, 
the hidden half-skyrmions align themselves, forming into 
a paired half-skyrmions line. With the mean field energy 
minimized, the deviation that we obtained is the conse- 
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quence of the cooperation of all the components of the 
condensate in the presence of the SO coupling and the ex- 
ternal trap. Thus, for spinor condensate with spin-orbit 
coupling strength n ~ 1 and comparable external traps, 
the usual plane wave phase or stripe phase are revised, 
bringing nontrivial detail structures. The present paper 
should help us to understand the cooperative roles of the 
SO coupling and the external potential in determining 
the ground state of the spinor condensates. 
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